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ABSTRACT 

The  application  of  numerical  optimization  techniques  to 
the  design  of  drive  shafts  is  demonstrated.   The  analysis 
investigated  the  effects  of  a  small  mass  imbalance  in  con- 
junction with  the  rotation  of  a  shaft  with  synchronous 
whirl. 
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I.   INTRODUCTION 

Optimization  and  composite  materials  are  two  fields 
that  have  seen  considerable  developments  in  recent  years. 
Advancements  in  these  technologies  have  provided  the  design 
engineer  with  new  dimensions,  the  depths  of  which  are  yet 
to  be  explored. 

The  major  material  characteristics  that  are  advanta- 
geously cultivated  in  composites  are  the  high  strength  and 
stiffness  to  weight  ratios.   The  weight  reduction  is  so 
significant  that  the  aircraft  and  aerospace  industries  have 
concentrated  a  major  portion  of  their  reserach  efforts  in 
this  field.   Investigations  in  improving  the  resistance  of 
composites  against  environmental  factors  such  as  temperatures, 
corrosion,  and  wear  are  continuing  and  findings  affirm  the 
contention  that  the  future  of  this  field  is  bright.   However, 
other  industries  are  slow  in  utilizing  these  materials  in 
basic  machinery  elements,  mainly  due  to  a  major  drawback — 
the  high  cost  of  manufacturing  and  processing.   It  is  in 
this  vein  that  the  importance  of  an  optimized  design  becomes 
apparent.   Advancements  in  fabrication  and  processing  of 
composite  materials  coupled  with  the  capability  to  optimize 
designs  may  well  lead  composites  into  a  more  competitive 
market. 
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with  the  ultimate  goal  of  optimizing  the  design  of  a 
drive  shaft  made  of  com.posite  materials,  a  research  program 
v/as  initiated.   This  report  is  the  first  in  this  investiga- 
tive effort. 

Before  the  intricacies  introduced  by  the  directional 
material  properties  of  composites  are  included  in  the 
analysis,  it  was  deemed  necessary  to  probe  the  characteris- 
tics of  the  design  of  an  isotropic  material.   This  will 
enable  one  to  foresee  the  trends  and  the  adaptive  measures 
necessary  when  the  composite  material  properties  are  con- 
sidered.  The  application  of  numerical  optimization  tech- 
niques to  the  design  of  an  isotropic  drive  shaft  is  demon- 
strated here.   The  effects  of  a  small  m^ass  imbalance  on  the 
bending  moment  and  deflection  of  a  shaft  with  synchronous 
whirl  are  investigated.   Mainly,  however,  the  objective  is 
to  develop  an  analysis  procedure  that  may  be  adapted  to  the 
peculiarities  of  composites  and  ultimately  be  used  in  op- 
timizing a  composite  drive  shaft  design. 

Section  2  defines  the  scope  and  lim.itations  of  the 
analysis  and  optimization. 

Section  3  discusses  the  analysis  itself,  the  equations 
and  the  modes  of  failure  that  were  considered.   Examples  of 
the  analysis  are  shown. 

Section  4  briefly  describes  the  optimization  program 
[1,2]  and  presents  the  results  of  optimizing  a  number  of 


shaft  designs.   Two  variable  function  space  diagrams  are 
provided  to  illustrate  the  design  space  for  different  load 
conditions. 

Section  5  contains  recommendations  for  future  investi- 
gations. 

Appendix  A  contains  a  description  of  the  FORTRAN  program 
used  in  the  analysis  and  the  program  text.  Appendix  B  shows 
the  derivation  of  equations  for  deflection,  moment  and  shear 
force. 
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II.   SCOPE  AND  LIMITATIONS 

The  analysis  applies  to  shafts  made  of  isotropic  ma- 
terials with  simply  supported  or  fixed-fixed  end  conditions. 
Basic  shaft  design  formulas  are  used  except  that  the  de- 
flection and  moment  equations  are  derived  to  include  the 
effects  of  a  small  mass  imbalance.  The  shaft  is  assumed 
to  rotate  with  synchronous  whirl. 

To-  insure  a  continuous  design  field,  the  optimization 
is  constrained  to  subcritical  speeds  only.   Supercritical 
speeds  introduce  feasible  design  regions  that  are  disjoint 
from  the  primary  design  space.   In  this  investigation,  these 
secondary  areas  are  considered  infeasible. 

Since  the  optimization  tends  to  a  large  radius-small 
thickness  design,  radial  stresses  are  neglected.   Hov/ever, 
buckling  of  thin  cylinders  due  to  torsion  and  compression 
is  considered.   The  equations  used  for  these  failure  criteria 
do  not  include  dynamic  effects.   A  constant  torque  and  a 
steady,  uniform  axial  force  are  the  loads  considered  in 
the  formulas  [3]. 

The  analysis  incorporates  the  capability  to  design  a 
shaft  that  may  be  used  in  tv/o  or  more  loading  conditions. 
For  example,  a  shaft  may  be  designed  to  transmit  50  HP  at 
300  RPM  with  an  axial  load  of  2000  lbs.  as  well  as  to  transmit 
150  HP  at  3000  RPM  with  a  1000  lb.  axial  load. 
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III.   ANALYSIS 

The  problem  is  a  shaft  required  to  transmit  a  specified 
horsepower  at  a  given  speed.   There  may  be  an  axial  load,  F 
or  an  internal  pres^efi^e/  P.  (See  Figure  1.)   Basic  equations 
are  used  throughout  the  analysis  and  fundamentals  such  as 
formulas  for  area,  volume,  moment  of  inertia,  polar  moment 
of  inertia  and  radius  of  gyration  are  not  repeated  here. 
However,  certain  stress  calculations  will  be  emphasized  to 
illustrate  the  flow  of  analysis. 

Appendix  B  shows  the  derivations  of  the  equations  for 
bending  moment,  deflection  and  shear  for  pinned-pinned  and 
clamped-clamped  end  conditions.   These  derivations  have 
taken  into  consideration  the  effects  of  a  small  mass  imbalance 
in  conjunction  with  the  rotation  of  a  shaft  with  synchronous 
whirl. 

The  inch-pound-second  system  of  units  (IPS)  is  used  in 
the  computer  program  analysis.   However,  by  changing  the 
equation  for  torque  (line  50  in  the  FORTRAN  text.  Appendix  A) 
to  reflect  the  International  System  of  Units,  the  program 
may  be  used  effectively  with  consistent  SI  units. 
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Figure  1.   Applied  Loads:   Axial  load  (F) ,  Internal  pressure  (P) , 

Torque  (T) ,  Mass  imbalance  (e) , 


TTX 


a.   Pinned-pinned  end  condition 


b.   Clamped-clamped  end  condition 


Figure  2.   Assumed  deflected  shape  of  the  shaft 
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A.   SIMPLY  SUPPORTED  BEAM 

The  deflected  shape  of  the  beam  is  assumed  as: 
6=6      sin  -T-  (see  Fig.  2) .   The  solution  of  the  govern- 
ing eqndLluniT  equations  is  given  in  Appendix  B.   The  deflection 
is  maximum  at  x  =  Z/2 


^max  =  384  ^i-k^  (1)  4  -'f  (1)  ^  '  '^•^^* 

2 

where;   K,  =  pAo) 

p  =  mass  per  unit  volume 

A  =  cross  sectional  area 

CO  =  shaft  speed  is  radians  per  second 

K  =  pAg  =  weight  per  unit  length 

e  =  eccentricity  of  mass  v;ith  respect  to  the  axis 
of  rotation 

i  =   shaft  length 

E  =  Young's  Modulus 

I  =  Moment  of  Inertia 

F  =  Axial  Load 
Certain  characteristics  of  this  equation  are  of  particular 
interest.   If  there  is  no  rotation  nor  axial  load,  the  formula 
reduces  to  the  classical  equation  for  maximum  beam  deflection. 

4 

^^^  =  384EI  ^^'^^^ 
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If  there  is  no  axial  load,  the  equation  becomes  unstable 

2  £  4 
(deflection  goes  to  infinity)  v/hen  EI  =  pAoo  (— )   or,  solving 

for  the  critical  speed. 


-   2  /  EI  ,-,  -,  N 

'^r        ^   \    4  (3.1c) 

^     V  pAil^ 


which  defines  the  first  fundamental  frequency  of  a  simply 

supported  shaft. 

£  2 
If  there  is  no  rotation,  instability  arises  when  EI  =  F (— ) 

IT 

or 

^c  "  ^^~¥'  ^^'^^^ 

I 

recognizable  as  Euler's  column  buckling  criteria. 

With  both  rotation  and  axial  load,  the  critical  speed 
becomes 


2  l^^ZJ^ 
to   =  7T   J (3.1e) 

To  avoid  the  instability  regions,  the  denominator  of 
equation  (3.1a)  must  be  given  a  lower  bound  greater  than  zero 
(point  A  in  Fig.  3) .   Additionally,  a  maximum  deflection  has 
to  be  imposed  to  prevent  computer  overflow  (point  B  in  Fig.  3) . 
These  limits  have  no  relevance  with  the  deflection  constraint 
used  in  the  design  optimization  process  (point  C  in  Fig.  3) , 
but  are  used  only  to  prevent  numerical  ill-conditioning  on 
the  computer. 
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^Deflection 


Denominator 


Figure  3.   Liinits  on  the  deflection.   The  solid  curve  repre- 
sents the  deflection  values  used  in  the  analysis. 
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The  bending  moment  corresponding  to  the  deflection  calcu- 
lated by  equation  (3.1a)  at  x  =  £/«  is 

»max  =  '^l  [|^  «(?)^'   +F6   +^  (3. If) 

When  the  weight  tends  to  cancel  the  deflection,  the  corre- 
sponding moment  at  the  same  cross  section  becomes 

2 

M  .   =  M     -. —  (3.1g) 

min    max     4  ^        ^i 

The  shear  force  is  maximum  at  x  =  0. 

K,6il   K,e£   K^Jl 

V^^^  =  -^^—-   +  -V-  +  -4-  (S.lh) 

max     TT      2      2 


B.       CLAMPED-CLAMPED    BEAr4 

The  deflected  shape  of  the  beam  is  assumed  as: 

6  =  -^  [  1  -  cos  ^^1  .   (See  Fig.  2b.)   The  deflection  is 


max 


imum  at  X  =  ^/ 'y- 


384  EI  -  0.746  K, £^  -  96  (-)   F 

1  7T 

As  in  the  simply  supported  beam,  similar  observations 
regarding  special  loading  conditions  are  apparent.   If  there 
is  neither  rotation  nor  axial  load,  the  formula  for  maximum 

beam  deflection  results 

4 
K  £ 

6    =  _,^  ^^  (3.2b) 

max    384  EI 
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with  no  axial  load,  instability  occurs  when 
384  EI  =  0.746  pAoo^  or 


0)^  =  22.68   J^^  (3.2c) 


v/hich  is  the  first  fundamental  frequency  for  a  fixed- fixed 

beam  based  on  the  assumed  deflection  used  here. 

When  there  is  no  rotation,  instability  arises  when 
2 
384  EI  =  96  ^-  or 

■^^ 

2 

F^  =  5 — ,  (3. 2d) 

C         n2. 


Euler's  column  buckling  criteria. 

With  rotation  and  axial  load,  the  instability  occurs  at 
a  smaller  speed  than  that  given  by  equation  (3.2c)/ 


oo  CO         \  E^  "  0.02533  F£^  ,.  ^. 

03^  =  22.68    / ;; (3.2e) 

^  ^'      pA£^ 


The  bending  moment  at  x  =  il/.  corresponding  to  the 
deflection  given  by  equation  (3.2a)  is 

2 

M  ^  -^      L  ^  +  ^  +   eil]   +  F6    +  -It-  (3.2f) 

max  24  2  7T  24 

A  lesser  bending  moment  occurs  at  the  same  cross  section 
when  the  weight  tends  to  cancel  the  deflection . 
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2 

M  .   =  M    =-Tr-  (3.2g) 

mm    max    12  v    ^/ 


The  maximum  shear  occurs  at  x  =  0. 


K^SZ        K   el        K^l 


The  preceding  calculations  are  performed  in  Subroutine 
BEND,  given  in  Appendix  A. 

C.   STRESS  ANALYSIS 

With  the  deflection,  shear  and  bending  moments  known  at 
the  critical  cross  sections,  the  analysis  proceeds  v/ith  the 
calculation  of  the  stresses  at  the  critical  stress  elements. 
At  midspan,  the  critical  stress  elements  are  on  the  outer 
surface  of  the  shaft  where  the  axial  stress  (tensile  or 
compressive)  is  maximum  (Points  A  and  B,  Fig.  4) .   Since  the 
shaft  is  in  synchronous  whirl,  the  same  stress  element.  A, 
will  generally  tend  to  be  in  tension  throughout  the  rotation, 
with  the  magnitude  of  the  axial  stress  varying  as  the  weight 
effect  becomes  alternatingly  additive  and  subtractive. 
Element  A  is  in  the  direction  of  the  mass  eccentricity  from 
the  center  of  the  shaft  while  element  B  is  on  the  outer 
surface  in  the  opposite  direction. 

Most  of  the  following  equations  are  found  in  basic  design 
and  applied  mechanics  books.   References  will  be  noted  when 
necessary. 
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Axial  stress  due  to  bending  moment  (tension/compression) 

M    r 
a^  (max)  = = (3.3a) 

M  .   r 
/  .  ^     mm   o 
a^(min)  = 

Axial  stress  due  to  axial  load  (compression) 

a   =  ^  (3.3b) 

X    A 

The  axial  stress  due  to  internal  pressure  for  thin-wall 
cylinders  is: 

Pr. 
^x  =  TF  ^^-^^^ 

Superposition  with  careful  regard  to  sign  conventions 
provides  the  maximum  and  minimum  tensile  and  compressive 
axial  stresses. 

Hoop  stress  is  primarily  due  to  internal  pressure  and 

rotational  inertia  [4] 

Pr 
0^=^  (3.3d) 

^"^   ^Q    3+v  [1  +  (i^)(-A)2]  (3.3e) 

a,.  =  .  —r—    "•      3+v  ''r    -• 

H     g      4  o 

where  Y  =  weight  per  unit  volume 

0)  =  rotation  in  radians  per  second 
V  =  Poisson's  ratio 

r.  =  inside  radius 

1 

r  =  outside  radius 
o 

g  =  acceleration  due  to  gravity 
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The  shear  stress  at  x  -  %/ ^    is  caused  only  by  torsion 
(shear  force  is  zero) . 

Tr 
T^  =  ^  (3.3f) 

where  t_,  =  shear  stress  due  to  torsion 
T  =  torque 

J  =  polar  moment  of  inertia 
At  X  =  0,  the  critical  element  is  where  the  shear  stress 
due  to  the  shear  force  along  the  cross  section  is  additive 
with  respect  to  the  shear  stress  due  to  torsion  (element  C, 

Fig.  3) .   The  shear  stress  due  to  V    for  thin  wall  cylinders 

max  -^ 

[5]  is 

2  V 

max  /  -)  -)  \ 

^v  =  — A—  ^^-^^^ 

Superposition  of  equations  (3.3f)  and  (3.3g)  gives  the  total 
shear  stress  at  the  critical  stress  element. 

The  axial  stress  at  element  C  due  to  bending  moment  is 
zero.   The  axial  stress  at  this  point  is,  therefore,  caused 
only  by  the  axial  load  and  the  internal  pressure,  equations 
(3.3b)  and  (3.3c). 

The  hoop  stress  is  calculated  by  equations  (3.3d)  and 
(3.3e) . 
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stress  element  A 


stress   eler.ent  3 
a)      At  midspan 


b)       At   X   =    0 


stress  element  G 


o„    (hoop   stress) 

ri 


T    (shear   stress) 


a      (axial   stress) 

X 


c)   Typical  stress  element, 
stresses  are  neglected, 

Figure  4.   Critical  Stress  Elements 


Radial 
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D.   MODES  OF  FAILURE 

Eight  modes  of  failure  are  considered  in  this  investiga- 
tion.  The  first  two  evaluate  the  static  strengths  at  the 
critical  stress  elements  using  Von  Mises  criteria  [5] .   At 
midspan. 


%id  =rn  -  ^H  ^x  -^  ^x^  -^  ^  4  (^-^^^ 


At  X  =  0 


-/ 


4    -   "e   "k   "■   "l   ^    ^    ^T+V  <3-^'^> 


Failure  is  assumed  if  the  yield  strength  of  the  material  is 
equal  to  or  less  than  these  Von  Mises  stresses  multiplied  by 
a  safety  factor. 

A  limit  on  the  deflection  is  imposed.  A  common  practice 
is  to  specify  a  maximum  deflection  in  inch  per  foot  of  shaft 
length.   For  example,  as  used  in  this  report, 

'^max  ^  0-005  inch/foot  length 

Since  the  optimization  may  lead  to  a  thin  cylinder  design, 
buckling  due  to  torsion  and  compression  are  considered.   The 
critical  torque  and  critical  compressive  stress  for  buckling 
of  thin  cylinders  [3]  are; 


^o  =    ^Vo   75  •  4^  "•^  =  * 
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^   =    ^^  (3.4d) 

^    r^3(l-v2) 

Experimental  data  indicate  that  actual  failure  of  the  cylinders 

usually  occurs  below  fifty  percent  of  the  values   calculated 

from  these  equations.   To  accommodate  this  discrepancy  and 

the  fact  that  the  dynamic  effects  are  not  considered  in  the 

derivation  of  these  equations,  comparatively  higher  safety 

factors  must  be  used  for  these  failure  modes. 

The  shaft  is  also  checked  for  column  buckling  using  Johnson's 

or  Euler's  equations  [5]  as  applicable.   The  actual  slenderness 

ratio  is  compared  with  the  slenderness  ratio  where  transition 

occurs  to  determine  which  formulas  apply.   The  calculated 

critical  load,  using  the  appropriate   equations  for  the  end 

condition,  is  compared  to  the  actual  load,  again  incorporating 

a  safety  factor.   Equations  (3. Id)  or  (3. 2d)  are  used  if  the 

slenderness  ratio   falls  in  the  Euler  range.   In  the  Johnson 

region,  the  following  equations  [5]  apply. 


^c  _  ^   ,  fi.2  {3.4e) 


where  S   =  yield  strength  in  compression 
{=r)    =  slenderness  ratio 

,    _    r   C-,  2    1 

^  -  %^  He 

E  =  Young ' s  Modulus 
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n  =  1  for  simple  supports 

4  for  fixed-fixed  ends 

The  effects  of  the  fluctuating  axial  stresses  are  evaluated 

for  fatigue  failure  using  the  Goodman  diagram  [5] .   Von  Mises 

stresses  for  the  mean  and  alternating  stresses  are  calculated 

as  follows: 

o"_  _   +  a_  .  _ 

(3.4f) 


(3.4g) 


a     = 
m 

(J 
max 

+ 

a    . 
mm 

2 

a 
max 

a    . 
mm 

2 

2  2        2 

a'  =,/a„  -   o„   a     +   o        +  3  t_ 
m   V  H     H   m    m       T 


(3.4h) 


(3.4i) 
a      a 

The  Goodman  diagram  is  drawn  (Fig.  5)  and  the  safety  factor 
is  determined  analytically. 


Figure  5.   Goodman  Diagram 
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SF  =  5^  =    ,   ^^ (3.4j) 

m  r   a    Se-| 

m 


where  Se  =  fully  corrected  endurance  limit 
Su  =  ultimate  strength 

This  safety  factor  is  required  to  be  greater  than  a  specified 
value. 

Lastly,  failure  is  assumed  if  a  specified  percentage  of 
the  critical  speed  as  calculated  by  equations  (3.1e)  or  (3.2e) 
is  less  than  the  shaft  speed. 

E.   ANALYSIS  EXAMPLES 

The  examples  are  chosen  to  illustrate  several  aspects  of 
the  shaft  design  problem.   The  first  is  a  design  which  vio- 
lates almost  all  of  the  constraints  while  the  second  example 
is  one  within  the  feasible  region  of  the  design  field. 

The  material  used  in  all  the  examples  in  this  report  is 
UN  G10150  HR  steel,  the  properties  of  which  are  listed  in 
Table  I  [5] . 
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TABLE  I 

MATERIAL  PROPERTIES 

UNS  G10150  HR  STEEL 

Young's  Modulus  30. OE  06  psi 

Shear  Modulus  11. 5E  06  psi 
Poisson's  Ratio  0.292        _ 

Weight  Per  Unit  Volume  0.282  lbs/in 

Yield  Strength  27. OE  03  psi 

Ultimate  Strength  50. OE  03  psi 
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1.   Analysis  Example  1 

A  simply  supported  shaft  with  0.1  inch  eccentricity 
is  to  transmit  150  HP  at  3000  RPM  (3150  in-lb  torque) , 
with  an  axial  load  of  2000  lbs.    For  a  shaft  length  of 
120  inches,  inside  radius  of  1.0  inch  and  thickness  of 
1.0  inch,  the  design  violates  the  static  and  dynamic 
strength,  deflection,  buckling  due  to  compression,  and 
the  speed  constraints.   The  analysis  of  this  design  is  given 
below  where  an  asterisk  (*)  denotes  violated  constraints  on 
the  design. 


Volume  0.11310E 

Weight  0.31893E 

Moment  of  Inertia  0.11781E 

Polar  Moment  of  Inertia  0.23562E 

*Critrical    Speed  0.14767E 

Critical  Axial  Load  0.18764E 

Critical  Buckling  Torque  0.47469E 

*Critical  Buckling  Stress  0.905493 

*Def lection  0.24000E 

Bending  Moment  0.36420E 

Shear  Force  0,62722E 

Axial  Stress  at  A  0.95779E 

Axial  Stress  at  B  -0.95783E 

Hoop  Stress  0.25013E 

Torsional  Shear  Stress  0.26738E 

Total  Shear  Stress  0.13337E 

Mean  Axial  Stress  0.95771E 

Alternating  Axial  Stress  0.81250E 

*Von  Mises  at  Midspan  0.95785E 

*Von  Mises  at  x  =  0  0.23100E 


04  in 

03  lbs 
02  in4 
02  in^ 

04  RPM 

06  lbs 
08  in  lb 

07  psi 

02  in 

08  in  lb 

06  lb 

07  psi 
07  psi 

03  psi 
03  psi 

06  psi 

07  psi 
03  psi 
07  psi 
06  psi 
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The  calculated  deflection  is  the  maximum  value 

imposed  to  avoid  numerical  ill-conditioning  (point  B,  Fig.  3) 

Since  the  deflection  affects  most  of  the  relevant  parameters 

in  the  analysis ,  the  calculated  stresses  are  unreliable  and 

not  valid.   These  values,  however,  tend  to  direct  the  design 

towards  the  feasible  region. 

2 .   Analysis  Example  2 

This  shaft  is  similar  to  the  first  example  in  design 

conditions  and  loading.   The  only  difference  is  that  the 

dimensions  are  such  that  no  constraints  are  violated,  but 

two  constraints  are  active  or  critical.   The  shaft  is 

120  inches  long  with  an  inside  radius  of  6.39  inches  and 

thickness  of  0.0369  inch.   This,  in  fact,  is  the  optimized 

design  for  the  given  loading  as  calculated  by  COPES-CONMIN 

[1,2].   The  results  of  the  analysis  follows,  where  a  double 

asterisk  (**)  denotes  values  at  their  design  limit. 

Volume  0.17866E  03  in^ 

Weight  0.50382E  02  lbs 

Moment  of  Inertia  0.30629E  02  in"^ 

Polar  Moment  of  Inertia  0.61258E  02  in"^ 

**Critical  Speed  0.60061E  04  RPM 

Critical  Axial  Load  0.39557E  05  lbs 

**Critical  Buckling  Torque  0.31490E  05  in  lb 

Critical  Buckling  Stress  0.10400E  06  psi 

Deflection  0.43813E-01  in 

Bending  Moment  0.17026E  05  in  lb 

Shear  Force  0.84895E  03  lbs 

Axial  Stress  at  A  0.22326E  04  psi 

Axial  Stress  at  B  -0.49193E  04  psi 

Hoop  Stress  0.29778E  04  psi 

Torsional  Shear  Stress  0.33079E  03  psi 

Total  Shear  Stress  0.14712E  04  psi 

Mean  Axial  Stress  0.20739E  04  psi 

Alternating  Axial  Stress  0.15872E  03  psi 

Von  Mises  at  Midspan  0.69313E  04  psi 

Von  Mises  at  x  =  0  0.39193E  04  psi 
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In  each  of  the  examples,  the  analysis  proceeds  as 
follows.   All  shaft  characteristics  like  cross  sectional 
area,  volume,  weight,  moment  of  inertia  and  polar  moment 
of  inertia  are  first  calculated.   Parameters  that  depend 
only  on  geometry  such  as  critical  torque  and  compressive 
stress  for  buckling  of  thin  cylinders,  and  critical  speed 
are  then  determined  before  Subroutine  BEND,  where  the  de- 
flection, bending  moments  and  shear  force  at  the  critical 
cross  sections  are  computed,  is  called. 

At  midspan,  the  absolute  values  of  the  axial 
stresses  at  elements  A  and  B  (see  Fig.  4a)  are  compared 
to  determine  which  element  is  critical.   In  both  examples, 
element  B  is  considered  critical  for  strength  evaluation. 
At  X  ss  0,  the  critical  element,  C,  is  where  the  torsional 
shear  stress  and  the  shear  stress  due  to  the  shear  force 
along  the  cross  section  are  additive.   The  Von  Mises 
stresses  for  each  of  the  critical  elements  are  then  calcu- 
lated using  equations  (3.4a)  and  (3.4b).   The  rest  of  the 
analysis  is  a  straight  forward  application  of  the  equations 
for  the  critical  parameters  and  the  constraints. 
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IV.   OPTIMIZATION,  RESULTS  AND  CONCLUSIONS 

A.   OPTIMIZATION  AND  RESULTS 

There  are  a  number  of  optimization  programs  available/ 
each  one  using  different  techniques  in  locating  the  desired 
optimiim  design.   COPES-CONMIN  [1/2]  is  a  versatile  program 
that  may  be  used  for  sensitivity  analysis  and  two  variable 
function  space  study  as  well  as  an  optimization  tool. 
Provided  with  a  user  supplied  analysis  program  (Subroutine 
ANALIZ)  where  the  objective  function/  constraints  and  other 
relevant  parameters  are  calculated,  it  determines  a  usable 
and  feasible  sector  from  which  a  search  direction  is  chosen. 
This  choice  of  search  direction  is  made  by  using  such  informa- 
tion as  the  gradients  of  the  objective  function,  gradients  of 
active  and  violated  constraints  and  the  degree  of  push-off 
these  constraints  violations  generate.   This  iterative  process 
of  minimizing/maximizing  the  objective  function  by  changing 
the  design  variables  is  terminated  when  no  further  improvement 
can  be  made. 

In  this  exercise,  the  objective  function  to  be  minimized 
is  the  volume,  with  the  radius  and  thickness  as  design 
variables.   Several  design  examples  are  shown  to  demonstrate 
the  effects  of  the  different  parameters  on  the  optimized 
design  and  on  the  design  space  itself.   Seven  examples  are 
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presented  corresponding  to  the  loading  conditions  as  listed 
in  Table  II.   Two  variable  function  space  diagrams  are 
generated  to  illustrate  these  effects.   The  numbered  curves 
in  the  graphs  represent  the  corresponding  constraints  shown 
in  Table  III.   The  shaded  side  of  a  curve  represents  the  un- 
feasible area  of  the  design  space  with  respect  to  the  cor- 
responding constraint.  The  dashed  curves  represent  contours 
of  the  objective  function,  the  numbers  indicating  the 
volume  in  cubic  inches. 

TABLE  II 
OPTIMIZATION  EXAMPLES 


No. 

Loads 

Horsepower 

RPM 

Pressure 

Axial  Load 

(psi) 

(lbs) 

1 

150 

300 

0 

0 

2 

150 

300 

1000 

0 

3 

150 

300 

0 

2000 

4 

150 

3000 

0 

0 

5 

150 

3000 

1000 

0 

6 

150 

3000 

0 

2000 

7* 

150 

300 

0 

0 

150 

3000 

0 

2000 

*Example  No.  7  is  a  design  for  two  load  conditions. 
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TABLE  III 
DESIGN  CONSTRAINTS 
Number  Constraint 

1  Static  Strength  at  Midspan 

2  Static  Strength  at  x  =  0 

3  Deflection 

4  Buckling  due  to  Torque 

5  Buckling  due  to  Compression 

6  Column  Buckling 

7  Dynamic  Strength  (Fatigue) 

8  Speed 
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EX^^PL5    1.     :       INITUL    05SISM 
SHAFT    DI^iENSIDNS    : 

RI    =    O.IOOOOE    01      TH    =    O.IOOOOE    01       SL    =    0.12000E 
OESIGM    CONDirnNS    : 

ECCEMTRieiTV    =    0,100005    00 
END    CDNDITIO^JS    :    PIMNED-PI NMED 
LOADS    : 

HORSEPOWER                RPM                     PRESSURE  AXIAL    LOAD 

0.15000E    03      0.30000E    03      0,0  0.0 

^MALYSIS    RESULTS 

T0R9UE(1I  0.31500E    05 

VOLUME  0.11313E    04 

WEIGHT  0.31893P    03 

MOMENT    OF    INERTIA  0.11731E    02 

POLAR    MOMENT    OF    INERTIA  0.23562E    02 

criti:al  sPEn  o.i^szbe  o^ 

CRITICAL    AXIAL    LOAD  0.1875^E    06 

CRITICAL    BUC<LING   tgrqjc  0.47<f59E    08 

CRITICAL    BUC<LIN1G    STRESS  0.905't9E    07 

DEFLECTION  0.2658^E-01 

BENDING    MOMEvIT  0.50477E    OA- 

SHEAR    FORCE  0,2071^E    03 

AXIAL    STRESS    AT    A  0.85692E    03 
AXUL    STRESS    ^T    B                                   -0.85692E    03 

HOOP    STRESS  0.25013E    01 

TORSIONAL    SHEAR    STRESS  0.26738E    04 

TOTAL    SHEAR    STRESS  0.27178E    04 

MEAN    AXIAL    STRESS  0.447&OE    02 

ALTERNATING    AXIAL    STRESS  0.81216E    03 

VON    MISES    AT    MIDSPAN  0.47095E    04 

VON    MISES    AT    X=0.  0,47073E    04 
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EX^^^PLE   1.     :       OPTIMUM    DESIGN 
SHAFT    DIMENSID^JS    : 

RI    =    0.25847=    01       TH    =    0.1il24£    00      SL    =    0. 12000E 
DESIGN    CONDITIONS    : 

ECCENTRICITY    =    O.IOOOOE    00 
END    CDNDITIO^S    :    P I NNED-PI NNED 
LOADS     : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15000E    03      0.30000E    03       0.0  0.0 

^^liVLYSIS    RESULTS 

TORQJEd)  0.3i503E    05 

VOLUME  0.2214SE    03 

WEIGHT  0.62451E    02 

MOMENT    OF    INERTIA  0.64353E    01 

POLAR    MOMENT    3F     INERTIA  0.12371E    02 

CRITICAL    SPEED  0.247S7E    04 

CRITICAL    AXIAL    LOAD  0,45137E    05 

CRITICAL    BUC<LING   TORQJE  0.31500E    06 

CRITICAL    BUC<LING   STRESS  0.74729E    06 

DEFLECTION  0.92753E-02 

BENDING    MOMEvlT  0.95479E    03 

SHEAR   FORCE  0.39681E    02 

AXI4L    STRESS    ^T    A  0.39998E    03 
AXIAL    STRESS    AT    B                                  -0.39998E    03 

HOOP    STRESS  0.51655E    01 

TORSIONAL    SHE^R    STRESS  0.65980E    04 

TOTAL    SHEAR    STRESS  0.66iflOE    04 

MEAN    AXIAL    STRESS  0.75507E    01 

ALTERNATING    AXIAL    STRESS  0.39243E    03 

VON    MISES    AT    MIDSPAN  0.11435E    05 

VON    MISES    AT    X=0.  0.11503E    05 
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EXAMPLE    2.    :       IMITIAL    OESIG^J 
SHAFT    DIMENSI3MS    : 

RI    =    0.10000=    01       TH    =    O.IOOOOE    01       SL    =    0.12000E 
DESIGN    CONOiri'!)NS    : 

ECCENJTRICITV    =    O.IOOOOE    00 
END    C3NDITI0^<S     :    PI  NNEO-PI NNED 
LOADS     : 
HORSEPOWER  RP'A  PRESSURE  AXIAL    LOAD 

0.15000E    03      0.30000E    03       O.IOOOOE    0^      0.0 

^MALYSIS    RESULTS 

TORQJE(l)  0.31500E    05 

VOLU^iE  0.11313E    0^ 

WEIGHT  0.31893E    03 

MOMEMT    OF    I>JERTIA  0.11731=    02 

POLAR    MOMENT    DP    INERTIA  0.23562E    02 

CRITICAL    SPE=D  0.1^828E    04 

CRITICAL    AXIAL    LOAD  0.18764E    06 

CRITICAL    BUC<LING   TORQUE  0.474S9E    08 

CRITICAL    BUC<LING   STRESS  0.90549E    07 

DEFLECTION  0.26584E-01 

BENDING    MOMEMT  0.50477E    04 

SHEAR   FORCE  0.2071^E    03 

AXIAL    STRESS    AT    A  0,18569E    04 

AXIAL    STRESS    AT    B  0.14303E    03 

HOOP    STRESS  0.66917E    03 

TORSIONAL    SHEAR    STRESS  0.26738E    04 

TOTAL    SHEAR   STRESS  0.27178E    04 

MEAN    AXIAL    STRESS  0.10448E    04 

ALTERNATING    AXIAL    STRESS  0.81216E    03 

VON    MISES    AT    MIOSPAN  0.49093E    04 

VON    MISES    AT    X=0.  0.4754&E    04 
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EXi^lPLE    2.     :       DOTnUM    OSSIGM 
SHAFT    DIMENSIONS    : 

RI    -    0.18625E    01       TH    =    0.25560E    00       5L    =    0.12000E 
DESIGN    CONDITIONS    : 

ECCEMTRICITV    =    O.IOOOOE    00 
END   CONDITIONS    :     PI  NNED-PI  NnjEO 
LOADS     : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15000E    03       0.30000E    03       O.IOOOOE    0*^      0,0 

ANALYSIS    RESULTS 

TORQJE(l)  0,31500E    05 

VOLUME  0,39958E    03 

WEIGHT  0.11268E    03 

MOMENT    OF    INERTIA  0.6657BE    01 

POLAR    MOMENT    OF    INERTIA  0.13316E    02 

CRITICAL    SPE=0  0.1875^E    0^ 

CRITICAL    AXIAL    LOAD  0.751'+^E    05 

CRITICAL    BUC<LING    TORQUE  0.23553E    07 

CRITICAL    BUC<LING    STRESS  0.22602E    07 

DEFLECTION  0.16358E-01 

BENDING    MOMEvJT  O.I7^76E    0^ 

SHEAR    FORCE  0.722^7E    02 

AXIAL    STRESS    AT    A  0,75709E    0^ 

AXIAL    STRESS    AT    B  0.64537E    04 

HOOP    STRESS  0.65487E    04 

TORSIONAL    SHEAR    STRESS  0.503'+4F    04 

TOTAL    SHEAR    STRESS  0.50773E    04 

MEAN    AXIAL    STRESS  0.70306E    04 

ALTERNATING    AXIAL    STRESS  0.54027E    03 

VON    MISES    AT    MIOSPAN  0.11254E    05 

VON    MISES    AT    X=0.  0,10965E    05 
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EX^^PLE    3.     :       I^IITIAL    OtSISM 
SHAFT    DIMENSIONS    : 

RI    =    0.10000=    01       TH    =    O.IOOOOE    01       SL    =    0, 12000E 
DESIGN    CONDirnNS    : 

ECCEMTRICITY    =    O.IOOOOE    DO 
END    C3NDITI0MS    :     PI NNEO-PI N^ED 
LOADS     : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15300E    03       0.30000E    03      3.0  0.20000E    04 

4VJALYSIS    RESULTS 

TORQUECll  0.31500E  05 

VOLUME  0.11310E  04 

WEIGHT  0.31893E  03 

MOMENT    Of    INERTIA  0.11731E  02 

POLAR    MOMENT    OF    INERTIA  0.23562=  02 

CRITICAL    SPEED  0.14767E  04 

CRITICAL    AXIAL    LOAD  0.18754E  06 

CRITICAL    BUC<LING   TORQUE  0.47469E  08 

CRITICAL    BUC<LING    STRESS  0.90549E  07 

DEFLECTION  0.26814E-01 

BENDING    MOMENT  0.54144=  04 

SHEAR    FORCE  0.2072DE  03 

AXIAL     STRESS    AT    A  0.70597E  03 

AXIAL    STRESS    \T    B  -0.11314E  04 

HOOP    STRESS  0.25013E  01 

TORSIONAL    SHEAR    STRESS  0.26733E  04 

TOTAL    SHEAR    STRESS  0.27178E  04 

MEAN    AXIAL    STRESS  -0.10519E  03 

ALTERNATING    AXIAL    STRESS  0.81216E  03 

VON    MISES    AT    MIDSPAN  0.47&77E  04 

VON    MISES    AT    X=0.  0.47073E  04 
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EXAMPLE   3.     :       OPTIMUM    DESIGN 
SHAFT    DIMENSIONS    : 

RI    =    0.25907=    01      TH    =    0.1I119E    00      SL    =    0. 12000E 
DESIGN   CONDITIONS    : 

ECCEMTRICIiy    =    O.IOOOOE    00 
END    :DNDITIOvIS    :    PI  N^JED-PI  NNED 
LOADS    : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15D00E    03       0.30000E    03      0.0  0.20000E    0^ 

^viALYSIS    RESULTS 

TOROUE(l)  0.31500E    05 

VOLUME  0.22186E    03 

WEIGHT  0.6256VE    02 

MOMENT    OF     INERTIA  0.6'V766E    01 

POLAR    MOMENT    DF    INERTIA  0.129535    02 

CRITICAL    SPEED  0.2^63S=    0^ 

CRITICAL    AXIAL    LOAD  0.45241E    05 

CRITICAL    BUCCLING    T3R0JE  0.31500E    06 

CRITICAL    BUC<LING    STRESS  0.74527E    06 

DEFLECTION  0.93752E-02 

BENDING    MOMcvJT  0.9317LE    03 

SHEAR    FORCE  0.39758E    02 

AXIAL    STRESS    \T    A  -0.67221E    03 

AXIAL    STRESS    AT    B  -0.1^913E    04 

HOOP    STRESS  0.51838E    01 

TORSIONAL    SHEAR    STRESS  0.65705E    04 

TOTAL    SHEAR    STRESS  0.66136E    04 

MEAN    AXIAL    STRESS  -0.i0637E    04 

ALTERNATING    AXIAL    STRESS  0.39151E    03 

VON    MISES    AT    MIDSPAN  0.11478E    05 

VON    MISES    AT    X=0.  0.11455E    05 
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EX\MPLE    ^.     :       IMITIAL    DESIGVI 
SHAFT    DIMENSI3VIS    : 

RI    =    0.10000=    01      TH    =    O.IOOOOE    01       SL    =    0.12000E 
DESIGN    C0N0iri3NS    : 

ECCEMTRICITy    =    0.10000=    00 
END    CDNOITIQVJS    :    PI  NNED-PI NMED 
LOADS    : 
HORSEPOWER  RPH  PRESSURE  AXIAL    LOAD 

0.15000E    03       0.30000E    0^      0.0  0.0 

^^JALYS!S    RESULTS 

TORQJE(l)  0.31500E  04 

VOLU^Ir  0.11310E  04 

WEIGHT  0.31893E  03 

MO^EMT    OF    IMERTIA  0.11781E  02 

POLAR    MOf^ENT    OF    INERTIA  0.23562E  02 

CRITI:aL    speed  0.14828E  04 

CRITICAL    AXIAL    LOAD  0.1876^E  06 

CRITICAL    3UC<uING    TORQUE  0.47'*59E  08 

CRITICAL    BUC<LING    STRESS  0.90549E  07 

DEFLECTION  0.24000E  02 

BENDING    MOMEVJT  0.23331E  08 

SHEAR    FORCE  0.62722E  06 

AXIAL    STRESS    ^T    A  0.4040SE  07 

AXIAL    STRESS    AT    B  -0.40406E  07 

HOOP    STRESS  0.25013E  03 

TORSIONAL    SHEAR    STRESS  0.26733E  03 

TOTAL    SHEAR    STRESS  0.13337E  06 

MEAN    AXIAL    STRESS  0.40397E  07 

ALTERNATING    AXIAL    STRESS  0.81200E  03 

VON    MISES    AT    MIDSPAN  0.40404E  07 

VON    MISES    AT    X=0.  0.23100E  06 
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EXX^PLE    4,     :       3PTIMUM    05SIGM 
SHAFT    DIMENSIDMS    : 

RI    s    0.63T98E    01       TH    =    0.36963E-01       5L    =    0.12000E 
DESIGN    CONDITIONS    : 

ECCEMTRICITr    -    O.IOOOOE    00 
END    C3NDITI0MS     :    PIN^JEO-PI  NNEO 
LOADS    : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15000E    03      0.30000E    04      0.0  0.0 

AVIALYSIS    RESULTS 

TORQUE(l)  0.31500E  04 

VOLUME  0.17331E  03 

WEIGHT  0,50283E  02 

MOMENT    OF    I\|E=^TIA  0.30415E  02 

POLAR    MOMENT    DF    INERTIA  0.60833E  02 

CRITI:aL    speed  0.60005E  04 

CRITICAL    AXIA:.    LOAD  0.3947S5  05 

CRITICAL    BUC<LING    TORQUE  0.31494E  05 

CRITICAL    BUC<LING    STRESS  0.10432E  06 

DEFLECTION  0.43922E-01 

BENDING    MOMENT  0.76214E  04 

SHEAR    FORCE  0.84773E  03 

AXIAL    STRESS    AT    A  O.i6079E  04 

AXIAL    STRESS    AT    B  -0.16079E  04 

HOOP    STRESS  0.29628E  04 

TORSIONAL    SHEAR    STRESS  0.33228E  03 

TOTAL    SHEAR    STRESS  0.14733F  04 

MEAN    AXIAL    STRESS  0.14488E  04 

ALTERNATING    AKIAL    STRESS  0.159I3E  03 

VON    MISES    AT    MIOSPAN  0.26327E  04 

VON    MISES    AT    X=0.  0.39103E  04 
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EX^^OLE    5.     :       niTUL    OESIGVJ 
SHAFT    DIMENSIDNS    : 

RI    =    O.IOOOOE    01      TH    =    O.IOOOOE    01       SL    =    0,12000E 
DESIGiM    CONDiriDNS    : 

ECCENJTRIGIT/    =    0.10000=    00 
END    CDNDITIO^S     :    PINNED-PI NMED 
LOADS     : 
HORSEPOWER  RPM  PRESSURE  AXTAL    LOAD 

0.15000E    03      0.30000E    0^      O.IOOOOE    OV      0.0 

AMALYSIS    RESULTS 

TORQUEd)  0.31500E    0^ 

VOLUME  0.11313E    04 

WEIGHT  0.31893E    03 

MOMENIT    OF    IMERTIA  0.11731^    02 

POLAR    MOMENT    3F    INERTIA  0.23552E    02 

CRITKAL    SPE=3  O.l-VBZSE    04 

CRITICAL    AXIAL    LOAD  0.1876^E    06 

CRITICAL    BUCKLING   TORQUE  0.474S9E    08 

CRITICAL    BUCCLING    STRESS  0.905^9E    07 

DEFLECTION  0.24000E    02 

BENDING    MOMENT  0.23801E    08 

SHEAR    F3RCE  0.62722E    06 

AXIAL    STRESS    AT    A  0.40416E    07 

AXIAL    STRESS    AT    B  -0.40396E    07 

HOOP    STRESS  0,91680E    03 

TORSIONAL    SHEAR    STRESS  0.26738E    03 

TOTAL    SHEAR    ST^^ESS  0.13337E    06 

MEAN    AXIAL    STRESS  0.40407E    07 

ALTERNATING    AXIAL    STRESS  0.81200E    03 

VON    MISES    AT    MIDSPAN  0.40411E    07 

VON    MISES    AT    X=0.  0.2310DE    06 
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EX\MPLE    5.     :       OPTIMUM    DESIGN 
SHAFT    DIMENSI3MS    : 

RI    =    0.608^^E    01       TH    =    3.61266E    00       SL    =    0.12000= 
DESIGN    CONDITIONS    : 

ECCENTRICITV    =    O.IOOOOE    00 
END    CDNDITIOnIS    :    P INNED-PI  NnJEO 
LOADS     : 
HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

O.15D00E    03      0.30000E    0^      O.IOOOOE    0^      0.0 

^^^LYSIS    RESU.TS 

TORO'JEdl  0.31500E    04 

VOLUME  0.29521E    04 

WEIGHT  0.83250E    03 

MOMENT    OF    INERTIA  0.50353E    03 

POLAR    MOMENT    DF    INERTIA  O.IOOTIE    04 

CRITICAL    SPEED  0.60003E    04 

CRITICAL    AXIAL    LOAD  0.65357E    06 

CRITICAL    BUCCLING    TOROJE  0.34401F    08 

CRITICAL    BUCKLING   STRESS  0.16557F    07 

DEFLECTION  0.43926E-01 

BENDING    MOMEvJT  0.12619E    06 

SHEAR   FORCE  0.14035E    05 

AXUL    STRESS    4T    A  0.11609E    05 

AXIAL    STRESS    AT    B  0.825Z5E    04 

HOOP    STRESS  0.12589E    05 

TORSIONAL    SHE4R    STRESS  0.209^8E    02 

TOTAL    SHEAR    STRESS  0.11620E    04 

MEAN    AXIAL    STRESS  0.11443E    05 

ALTERNATING    AXIAL    STRESS  0.16609E    03 

VON    MISES    AT    ^IDSP^VN  0.12129E    05 

VON    MISES    AT    X=0.  0.12749E    05 
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EX\MPLE   6,    :       INITIAL    DESIGM 
SHAFT    DIMENSn^lS    : 

RI    =    0.10000:    01       TH    =    0. lOOOOE    01       SL    '    0.12000E 
DESIGN    CONDITIONS    : 

ECCENTRIGITy    =    O.IOOOOE    30 
END    CDNDITIOMS    :    P!  NNEO-PI NNEO 
LOADS     : 

HORSEPOWER                 RPM                      PRESSURE  AXIAL    LOAD 

0.15000E    03      0.30000E    04      0.0  0.20000E    04 

A^ALYSI S    RESULTS 

TORQJE(l)  0.31500E    04 

VOLUME  O.11310E    04 

WEIGHT  0.31893E    03 

MO^iEMT    OF    IVIE^TIA  0.11731E    02 

POLAR    MOMENT    DF    INERTIA  0.23552E    02 

CRITICAL    SPEED  0.14757E    04 

CRITICAL    AXIA.    LOAD  0.13764E    06 

CRITICAL    BUC<LING    TORQUE  0.47469E    08 

CRITICAL    BUC<LING   STRESS  0.905^9E    07 

DEFLECTIOM  0.2400DE    02 

BENDING    MOMEvIT  0.56420E    08 

SHEAR    FORCE  0.62722E    06 

AXIAL    STRESS    AT    A  0.95779E    07 
AXIAL    STRESS    AT    B                                   -0.95733E    07 

HOOP    STRESS  0.25013E    03 

TORSIONAL    SHEAR    STRESS  0.26733E    03 

TOTAL    SHEAR    STRESS  0.13337E    06 

MEAN    AXIAL    STRESS  0.95771E    07 

ALTERNATING   AXIAL    STRESS  0.81250E    03 

VON    MISES    AT    MIOSPAN  0.95735E    07 

VON    MISES    AT    X=0.  0.23100E    06 
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EX\MPLE    6.    :       JPTI^UM    DESIGVJ 
SHAFT    DIMENSIONS    : 

RI     =    0.63959E    01       TH    =    0.369^2E-01       SL  =    0.12000E 
DESIGN    CONOiriONS    : 

ECCENTRICITY    =    O.IOOOOE    DO 
END    C3N0ITI0NS    :     PINNED-PI NNEO 
LOADS     : 

HORSEPOWER  RPM  PRESSURE  AXIAL    LOAD 

0.15DOOE    03      0,30000E    0^      0,0  0.20000E    04 

AMALYSIS    RESULTS 

TORQJECl)  0.31500E  04 

VOLUME  0.i7866E  03 

WEIGHT  0.50382E  02 

MOMENT    OF    INERTIA  0.30629E  02 

POLAR    MOMENT    3F     INERTIA  0.61253E  02 

CRITICAL    SPEED  0.60051E  04 

CRITICAL    AXIAL    LOAD  0.39557E  05 

CRITICAL    BUC<L!NG    TOROUE  0,31490E  05 

CRITICAL    BUC<LING    STRESS  0.10400E  06 

DEFLECTION  0.438I3E-01 

SENDING    MOME^IT  0.17025E  05 

SHEAR    FORCE  0.84895E  03 

AXIAL     STRESS    AT    A  0.2232SE  04 

AXIAL     STRESS    A"^    B  -0.49193E  04 

HOOP    STRESS  0.29773E  04 

TORSIONAL    SHEAR    STRESS  0.33079E  03 

TOTAL    SHEAR    STRESS  0.14712E  04 

MEAN    AXIAL    STRESS  0.20739E  04 

ALTERNATING    AXIAL    STRESS  0.15372E  03 

VON    MISES    AT    MIDSPAN  0.69313E  04 

VON    MISES    AT    X=0.  0.39193E  04 
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EX^'^PLE    7.    :       niTUL    OESIGVI 
SHAFT    DIMENSIONS    : 

RI    =    O.IOOOOE    01       TH    =    O.IOOOOE    01       SL    =    0.12000E 
DESIGN    CONDiriGMS    : 

ECCENTRIGITy    =    0,10000=    30 
END    CDNDITia^S    :    PINNED-PI NNED 
LOADS     : 

HORSEPOWER                 RPM  PRESSURE              AXUL    LOAD 

0.L5300E    03       0.30000E    03  0,0                           0.0 

0.15000E    03       0,30000E    0^  0.0                           0.20000E    04 

AvlALYSIS    RESULTS 

TOROUE(l)  0.31500E    05 

T0R0JE(2)  0.31500E    OA- 

VOLUME  0.1131DE    04 

WEIGHT  0.31893E    03 

MOMENT    3F    INERTIA  0.11781c    02 

POLAR    MOMENT    OF    INERTIA  0.23562E    02 

CRITICAL    SPE=D  0.14767E    04 

CRITICAL    AXIAL    LOAD  0.18764E    06 

CRITICAL    BUC<LING    TORQUE  0.47469E    08 

CRITICAL    BUC<LING    STRESS  0-90549E    07 

DEFLECTION  0.24000E    02 

BENDING    MOMENT  0.56420E    08 

SHEAR    FDRCE  0.62722=    06 

AXIAL    STRESS    AT    A  0.95779E    07 

AXIAL    STRESS    AT    B  -0.95783E    07 

HOOP    STRESS  0.25013E    03 

TORSIONAL    SHEAR    STRESS  0.26733E    03 

TOTAL     SHEAR    STRESS  0.13337E    06 

MEAN    AXIAL    STRESS  0.95771E    07 

ALTERNATING    AXIAL    STRESS  0.81250E    03 

VON    MISES    AT    MIOSPAN  0.95785E    07 

VON    MISES    AT    X=0.  0.23100E    06 
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EXX^OLE    7.    :       DPTI^UM    OESIS^J 
SHAFT    0IMENSI3NS    : 

RI    =    0.63551E    01       TH    =    0.93009E-01       SL    =    0. 12000E 
DESIGN    CONOiriDiNJS    : 

ECCEMTRICITy    =    0.10000=    00 
END    C3NDITI0VIS     :    P INNED-PI NMED 
LOADS    : 

HORSEPOWER                 RPM  PRESSURE               AXIAL    LOAD 

0.15000E    03      0,30000E    03  0,0                           0.0 

0.15000E    03       0,30000E    0^  0.0                           0.20000E    04 

\M!\LYSIS    RESULTS 

TORQUE(l)  0.31500E    05 

T0RQJE(2)  0,31500E    04 

VOLUME  0.44893E    03 

WEIGHT  0.12660E    03 

MOMEMT    OF    nE^TIA  0.76650E    02 

POLAR    MOMENT    DF     INERTIA  0.15332F    03 

CRITICAL    SPEED  0.60000E    04 

CRITICAL    AXIAL    LOAD  0.99390E    05 

CRITICAL    BUC<LING    TORQUE  0.31571E    06 

CRITICAL    BUCKLING    STRESS  0.26122E    06 

DEFLECTIOM  0.43932E-01 

3EN0ING    MOMEvIT  0.42893F    05 

SHEAR    FORCE  0.21344E    04 

AXIAL    STRESS    AT    A  0,30733E    04 

AXIAL    STRESS    AT    B  -0.41425E    04 

HOOP    STRESS  0,29828E    04 

TORSIONAL    SHEAR    STRESS  0.132^BE    03 

TOTAL    SHEAR    STRESS  0.12736E    04 

MEAN    AXIAL    STRESS  0.29136E    04 

ALTERNATING    AXIAL    STRESS  0.15973E    03 

VON    MISES    AT    MIDSP^^N  0.62021E    04 

VON    MISES    AT    X=0.  0.37093E    04 
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B.   OBSERVATIONS  AND  CONCLUSIONS 

With  all  other  parameters  held  constant,  an  increase  in 
speed  drives  the  design  to  a  very  large  inside  radius  while 
reducing  the  thickness  considerably.   This  trend  is  so 
strong  that  a  reasonable  lower  bound  on  the  thickness  must 
be  given  to  prevent  numerical  ill-conditioning.   The  intro- 
duction of  eccentricity  increases  the  effect  of  rotation  on 
the  stresses  and  the  optimization  produces  slightly  larger 
values  for  the  design  variables.   This  effect  is  reduced, 
however,  by  the  constraint  on  the  deflection. 

As  the  internal  pressure  is  increased,  the  tendency  is 
to  produce  an  optimized  radius  that  is  smaller,  and  a  larger 
thickness.   This  may  be  explained  by  the  direct  effect  of 
the  internal  pressure  on  the  hoop  stress  which  in  turn  would 
require  a  larger  thickness  to  support  the  higher  stress. 

The  effect  of  axial  load  is  reflected  on  the  deflection. 
Again,  if  very  small  deflections  are  allowed,  the  axial  load 
effects  are  reduced.   The  major  parameter  that  is  affected 
is  the  axial  stress. 

A  notable  characteristic  shown  by  the  two  variable 
function  space  diagrams  is  the  tendency  of  the  strength 
constraints  to  curl  back  as  the  radius  is  increased.   This 
signifies  the  existence  of  an  upper  limit  on  size  for  a  given 
design  load.   The  increase  in  mass,  especially  at  high  speeds, 
increases  the  stresses  considerably  leading  to  violation  of 
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the  strength  constraints.  This  phenomenon  is  primarily  a 
result  of  the  eccentricity  of  the  mass  but  also  is  due  to 
the  static  weight  of  the  shaft  itself. 

It  is  also  noted  that  these  strength  constraints  generally 
run  close  to  each  other  in  an  almost  parallel  manner.   The 
dynamic  or  fatigue  constraint  maintains  a  position  such  that 
coming  from  the  feasible  range  towards  the  infeasible  area, 
it  is  activated  before  the  static  constraints  are  encountered. 
In  effect,  the  latter  constraints  are  redundant  in  rotating 
shafts.   However,  these  are  kept  in  the  analysis  to  allow 
solution  of  static  problems. 

The  versatility  of  the  program  should  not  be  lost  in  the 
emphasis  on  drive  shaft  applications.   With  the  appropriate 
loadings,  one  may  use  it  to  solve  structural  problems  like 
columns,  beams  and  pressure  vessel  designs.   The  objective 
function  need  not  be  confined  to  the  volume.   For  example, 
given  a  specific  initial  dimensions,  one  may  minimize  the 
Young's  modulus  to  determine  whether  or  not  a  certain 
material  may  be  utilized;  or  to  determine  how  low  the  modulus 
could  get  before  the  shaft  fails.   If  the  outer  dimension 
is  critical  to  the  design,  one  may  set  a  specific  value  on 
the  diameter  and  let  all  other  parameters  vary  during 
optimization . 

In  summary,  numerical  optimization  provides  an  efficient 
and  effective  way  of  improving  shaft  designs.   The  analysis 
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used  in  this  exercise  is  a  first  step  when  compared  to  the 
boundless  refinements  in  techniques  that  are  possible. 
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V.   RECOMMENDATIONS  FOR  FUTURE  ITWESTI CATIONS 

The  study  has  shown  the  feasibility  of  using  numerical 
optimization  techniques  in  the  design  of  drive  shafts  within 
the  limitations  imposed  in  the  analysis.   Further  studies  on 
the  same  design  field  may  be  pursued  by  eliminating  some  of 
these  limitations.   For  example,  the  unfeasible  region  defined 
by  the  speed  constraint  extends  only  to  some  discrete  distance 
and  a  supercritical  feasible  design  area  does  exist.   It  will 
be  interesting  to  develop  optimization  techniques  to  accommo- 
date this  disjoint  design  field  and  to  investigate  the  pecu- 
liarities this  will  introduce  to  the  analysis. 

Another  investigation  may  be  concerned  with  shafts  with 
variable  thickness  and  radius  along  its  length.   The  inclusion 
of  radial  stresses  in  the  calculations,  the  use  of  finite 
element  methods  in  three  dimensional  stress  analysis  of  the 
problem,  and  vibration  considerations  are  but  some  of  the 
aspects  of  research  that  may  be  useful  in  future  analyses  of 
composite  materials. 

The  challenge  of  including  the  unique  material  character- 
istics of  composites  in  the  design  will  invariably  stimulate 
a  lot  of  young  minds.   However,  the  handling  of  these 
direction  oriented  properties  in  stress  analysis  have  already 
been  investigated  [6  through  12],  and  the  tools  required  to 
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carry  out  the  program  are  available.   COMAND  [13]  offers  a 
simplified  composite  analysis  procedure  for  panels.   Such 
procedures,  adapted  with  the  program  used  in  this  investi- 
gation could  result  in  a  simplified  approach  to  composite 
drive  shaft  analysis.   In  this  procedure,  failure  of  the 
shaft  is  assumed  if  a  single  ply  fails  and  the  entire 
analysis  is  performed  in  a  ply  by  ply  basis.   Since  the 
stress  elements  have  to  be  in  a  differential  level,  the 
composite  material  is  assumed  to  be  homogenous  and  uni- 
directionally  isotropic  at  the  differential  range.   It  will 
also  be  necessary  to  assume  that  stresses  throughout  the 
thickness  of  the  ply  is  uniform  and  equal  to  the  stresses 
at  the  outer  element  of  the  ply.   Analysis  may  then  be 
performed  in  the  classical  procedure,  solving  for  the 
stresses  along  the  major  shaft  axes.   The  material  properties 
to  be  used  in  the  calculation  must  be  equivalent  values 
derived  from  the  directional  properties  using  transformation 
of  axes  equations.   References  10  and  13  present  procedures 
to  calculate  the  corresponding  material  properties,  stresses 
and  strains  along  the  major  axes  of  the  shaft  and  the  fiber 
orientation  coordinates. 

More  involved  investigations  could  deal  with  failure 
modes  of  the  whole  shaft  instead  of  a  single  ply,  three 
dimensional  analysis  of  plys  including  transverse  shear 
stresses,  delamination  problems,  thermal  effects,  and  others. 
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The  field  of  composites  also  opens  up  another  group  of 
considerations  like  the  effects  of  stacking  different  ply 
orientations  at  different  sequences,  unsyiranetrical  laminates, 
unequal  thicknesses  of  plys  and  other  factors  relevant  only 
to  composites. 
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APPENDIX  A 

The  FORTRAN  program  used  in  this  investigation  is 
described  here.   Subroutine  ANALIZ  provides  the  basic 
analysis  used  in  the  optimization.   It  reads  the  initial 
design  description  and  calculates  the  values  of  the  objec- 
tive function,  constraints,  and  all  other  parameters 
necessary  to  solve  the  analysis  problem.   COPES-CONMIN 
updates  the  design  to  minimize/maximize  the  objective 
function,  iterating  until  no  further  improvement  in  the 
objective  function  is  possible  without  violating  one  of 
the  constraints.   This  process  is  shown  in  the  flow  diagrams, 

The  global  catalog  lists  the  location,  FORTRAN  name, 
mathematical  symbol  and  description  of  the  parameters  used 
in  the  optimization.   These  parameters  are  contained  in  the 
labeled  COMMOM  block,  GLOBCM. 

There  are  at  least  five  input  cards  to  describe  the 
initial  design.   The  first  card  contains  the  initial 
dimensions — inside  radius,  thickness  and  shaft  length.   The 
second  card  describes  the  material  used.   Young's  modulus, 
shear  modulus,  Poisson's  ratio,  weight  per  unit  volume, 
yield  strength,  ultimate  strength  and  compressive  strength 
are  read  here.   The  endurance  limit  is  internally  calculated 
using  the  endurance  limit  factors  read  from  the  third  card. 
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The  size  factor,  CB,  is  internally  updated  as  the  shaft 
diameter  is  varied  during  optimization. 

The  fourth  and  fifth  cards  specify  the  design  and  load 
conditions.   The  parameter,  LC,  allows  the  user  to  choose 
the  end  conditions  applicable  to  his  requirements.   Pinned- 
pinned  and  clamped-c lamped  analyses  are  available  while 
clamped-pinned  and  cantilever  end  conditions  are  yet  to  be 
developed.   The  parameter  LCC  allows  multiple  loading  design. 
One  can  thereby  design  a  shaft  that  is  capable  of  supporting 
LCC  different  load  combinations. 

The  magnitude  of  mass  imbalance  may  be  read  as  a  specified 
distance  from  the  center  of  the  shaft,  EO,  or  as  a  specified 
fraction,  EC,  of  the  shaft  radius.   The  safety  factor 
generally  used  in  the  analysis  is  also  specified  in  the  fourth 
card. 

The  fifth  card  specifies  the  horsepower,  shaft  speed  in 
RPM,  internal  pressure  and  axial  load  for  each  loading. 
Additional  cards  are  used  up  to  LCC  number  of  loadings.   The 
torque  is  calculated  internally  in  inch-pounds.   The  follow- 
ing table  summarizes  the  required  information  for  subroutine 
ANALIZ. 
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TABLE  IV 
•INPUT  CARDS 
CARD  INFORMATION  FORMAT 

1  RI,  TH,  SC  8F10.2 

2  E,  G,  PR,  SPWT,  I/S 

SU,  YC  8E10.2 

3  CA,  CE,  CC,  CD,  CE, 

CF  8F10.2 

4  LC,  LCC,  EO,  EC,  SF       215,  3F10.2 

5  HP,  RPM,  P,  F  8F10.2 

Subroutine  BEND  calculates  the  bending  moment,  deflection 
and  shear  force  at  the  critical  cross-sections  of  the  beam. 
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SUBROUTINE  ANALIZ  (ICALC) 
FLOW  DIAGRAM 


ICALC  =  1 


READ   INITIAL  DIMENSIONS 
READ   MATERIAL  PROPERTIES 
READ   DESIGN  CONDITIONS 
PRINT  ALL  INPUT 

ICALC  =  2 

CALCULATE  GEOMETRIC  PARAMETERS 

CALCULATE  OBJECTIVE  FUNCTION 

CALL  SUBROUTINE  BEND 

.  CALCULATE  MAXIMUM  DEFLECTION 

.  CALCULATE  BENDING  MOMENTS 

.  CALCULATE  SHEAR  FORCE 

CALCULATE  STRESSES  AT  CRITICAL  STRESS 
ELEMENTS 

CALCULATE  CONSTRAINTS 
ICALC  =  3 


PRINT  ANALYSIS  RESULTS 
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OPTIMIZATION  FLOW  DIAGRAM 


INPUT 

1 

SUBROUTINE 
ANAL I Z 

SUBROUTINE 
BEND 

t 

r 

COPES 

^~ 

CONMIN 

i 

OUTPUT 
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TABLE  V 
GLOBAL  CATALOG  OF  PARAMETERS 


Global         FORTR?^N        Math 
Location Name   Symbol Definition 

Inside  Radius 

Thickness 

Shaft  Length 

Volume 

Torque 

Internal  Pressure 

Axial  Load 

Shaft  Speed 

Young's  Modulus 

Shear  Modulus 

Poisson's  Ratio 

Yield  Strength 
(Tension) 

Yield  Strength 
(Compression) 

Ultimate  Strength 

Outside  Radius 

Weight  Per  Unit 
Volume 

Moment  of  Inertia 

Polar  Moment  of 
Inertia 

Hoop  Stress 

Axial  Stress  at  A 

Minimum  Axial  Stress 

Axial  Stress  at  B 

Alternating  Axial 
Stress 

Mean  Axial  Stress 

Torsional  Shear 
Stress 


1 

Rl 

r . 
1 

2 

TH 

t 

3 

SL 

£ 

4 

VOL 

V 

5- 

•14 

TQ(IO) 

T 

15- 

■24 

P(10) 

P 

25- 

■34 

F(10) 

F 

35- 

■44 

RDS(IO) 

OJ 

45 

E 

E 

46 

G 

G 

47 

PR 

V 

48 

YS 

sy 

49 

YC 

s 

c 

50 

SU 

S 
u 

51 

RO 

r 
o 

52 

SPWT 

Y 

53 

AI 

I 

54 

AJ 

J 

55 

HS 

^H 

56 

AST 

^XA 

57 

ASTN 

X  (min) 

58- 

■67 

ASC(IO) 

^XB 

68- 

-77 

SALT (10) 

a    a  ' 
a,   a 

78 

SMEAN 

a 
m 

79 

SST 

^T 
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80 

81 
82 

83 


ssv 

SSMX 
TCRT 

SCSX 


max 


Shear  Stress  due  to 
Sheer  Force 

Total  Shear  Stress 

Critical  Torque  for 
Buckling 

Critical  Stress  in 
Buckling 


84- 

-93 

VMSI(IO) 

^m 

Von  Mises  (mean) 

94- 

-103 

VMS2 (10) 

^o' 

Von  Mises  at  x  = 

=  0 

104- 

-113 

VMS3(10) 

%id' 

Von  Mises  at  Midspan 

114 

WM 

'"c 

Critical  Speed 

115 

BM 

M 
max 

Maximum  Bending 

Moment 

116 

BMN 

M  . 
mm 

Minimum  Bending 

Moment 

117 

V 

V 
max 

Shear  Force 

118 

FC 

^c 

Critical  Column 
Load 

Buckling 

119 

SF 

SF 

Safety  Factor 

120- 

-129 

DEL(IO) 

6 

Deflection 

130- 

-239 

GC(10,10) 

GC 

Constraints 
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APPENDIX  B 


The  derivation  of  the  equations  for  the  deflection, 
bending  moments  and  shear  force  at  the  critical  cross- 


sections  of  the  shaft  are  presented  here.   Pinned[-pinned 
and  fixed-fixed  end  conditions  are  analyzed. 
1.   Pinned-Pinned 


-fc     X 


TTX 


Let   y   =   d      sin     „ 
-^  o  a 


(1.1) 


dQ   =   dM    (y+e)    co 

2 

dQ   =    pAoj       (y+e)    dx 

2 


Let   K^    =    pAw 


TTX 


dQ  =   K^    [d^   sin   -^  +   e]    dx 


(1.2) 


(1.3) 
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with  the  effect  of  weight: 

dQ  =  [K^(d^  sin  ^  +  e)  -  K^]  dx  (1.4) 


where  K2  =  pAg  =  numerically  equal  to  weight  per  unit  length 
*i  ^  d.£ 

Moment  at  x  =  x 


/*£  d  £ 

Q=   \   dQ  =  Kj_[^  _£_+e£]  -   K^    Z  (1.5) 


/ 


^xo  =  -  I  ^^-  ^o^  ^   ^  (^  -  ^o^  ^Q  (1-^) 


\  =  "  i  ""  ^1  ^^  ^"F  ^  1^  "■  t|  "  ^1^  ^"T  "■  ^^^'^   ^^'^^ 


1    2    ^.   ,  ,£.  2      7TX 
+  __  X  -  K3_  d^(-)   sm  — 


K2    2  ^  ^2^ 
-^   X   +  -^  X 


This  represents  the  moment  at  any  value  of  x.   If  the  axial 
load,  F,  is  included,  the  total  bending  moment  is 

M       =  M   +  IVL,  =  Ely" 
TOTAL     X    T     ^ 

where  M^  =  Fy  =  F  (d^  sin  ^)  (1.8) 

Integrating  twice: 

:iy  =  -4ir    [i^  X  -  2£x-^  +  x^]  +  K^  d   (-)  ^  sin  •—   (1.9) 

■^      24  1   O   7T  X, 


e; 


+  F  d   (-)^  sin  ^  -  ^  [i^x   -   2ix^   +   x^ ] 
O   TT         Jl     24  "■ 
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At  X  =   1/ ^  ^   y  =  d 

il    4  £    2         ^    ^1^^^         ^    S    ^ 

Eld      -   Kt    d      (-)^    -   Fd    (-)^   =   — r^T o4t- 

o  1      o      IT  o    IT  384  384 

^,4                      K   e    -   K^ 
^o   =   Im      t £    4      TT^  (^-^^^ 

1     77  77 


This  is  the  deflection  when  the  weight  effect  is  subtractive. 

The  maximum  deflection  occurs  when  the  weight  tends  to  increase 

the  deflection. 

5^4        K^e  +  K^ 
^max  "  384  ^ Tl Tl)  (1.11) 

1  TT  77 

Substituting  (1.10)  in  (1.7)  and  (1.8),  the  bending  moment 
at  X  =  £/«  is 

M  .   =  K,  [|  +  d  {-)     \      +  Fd ^  (1.12) 

mm    1  "-S    o   77'  ■■  8 

Again  the  maximum  moment  occurs  when  the  weight  effect  is 
additive 

p  p  2  K^£2 

M    =  K,  [|  +  d     {-)    \      +  Fd  +  -V-  (1.13) 

max    1  "-8    max   tt'  ■■  8 

From  Equation  (1.5),  the  maximum  shear  force  V    is  equal 
to  J  or 

K.  d    I         K,e£    K„Jl 

V     =  -1 ^^^      +  -i +  -1-  (1  14) 

^MAX       77      ^22  U.14; 
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2 .   Clamped -Clamped 

A  similar  derivation  procedure  for  clamped-clamped  ends 

using  a  shape  function: 

d  , 

y  =  -^  [1  -  cos  ^]  (2.1) 


and 


d    d 
dQ  =  K^  [-2.  -  ^   cos  ^  +  e]  dx  (2.2) 


Q  =  ^  dQ  =  K,il  [-^  +  e]  (2.3) 

o 


d 
c 


Moment  at  x  =  x 

o 


/»x 


M  =  M   -  §  X  +   1    (x   -  x)  dQ 
O    2       ^J       O 

o 


,  £   d  p^o   X  d     X  d       ^^ 

l/0,x,Tr    \    rOo     oo      2Trx 


=  M   - 

xd       xd       ^ 
+  x^e 2^  +   -^  cos  ^  -  xe]  dx      (2.4) 

M  at  any  value  of  x: 

K,  £   d  K,   d        ^ 

Ely"  =  M  =  M J-    {-J-  +   e)    x   +  -^    {-j-  +   e)    x^ 

K,d  £^  o 

-  2~  (1  ~  ^°^  ^~^  (2.5) 

87T 
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Integrating  once: 


K,       d  3  2,  K,d   l^x 


Ely'    =   M^   X   -   -^    (-^  +   e)     (V  -   ^    ^^  ^    "^ 


Stt 

1     O  „.  27TX      ,       _ 


Y*    =    0    at   X   =    0,    £/2    and  i,    c     =    0 

Stt 
Integrating  once  more: 

T7T  d  4  3,  2,2 

Kj_  ^  2  ^^     ^24  12  24    ^ 

d    ^"^  ^ 

o  2Trx    ,     -, 

-   -TTA      ^°^   -X-  ^   ^2 

32tt 

d    Z^ 
y=Oatx   =    0,    C2=  -^— ^ 

32tt 

^^  d  2  ^         d    l^  ^ 

Ely  /    o    ,       .      ,x    >   ,  ,v2     ,       o  ,,  27rx.        , r.    >.x 

^T^  =    (-2"  +   e)     (24-)  (x  -  £)       +   ^    (1    -    cos   -J-)        (2.6) 

1  32tt 


Y   =    d      at   X   =    £/„ : 
■^  o  2 


K^      o  2  96  4  3_g^4 


Eld           d    ^"^  o4  d    £^ 
o  _      o  e£  o 

K,  ~     768  384  ,^    4 

1  167T 
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4 
%    =  i-^ (2.7) 

384EI  -^4—  (^  +  24) 

With  axial  load  F: 

Fd  ^ 

Moment  due  toF  =  M^  =  M   -F    =M   -  -^    (1  -  cos  =^) 

F     o     Yp     o     2  i 

Integrating  once:  (2.8) 

Ely/  =  „_^  K  -  ^  (X  -  ^  Sin  2p)  +  Ci 


Yp'  =  0  at  X  =  0,  C,  =  0 

Fd 
Yp  =  0  at  X  =  e ,  M^  =  -^ 

where  M  =  Moment  at  the  ends  due  to  the  axial  load 
o 

Fd  2  Fd  2  .2  ^ 

tr-r  o      ,X     .  O     rX  ,        ^  2TrxT       ,      ^ 

EiY  =  -2-   (-T^    -  ^-   tx     ^  -T  ^°^  -T-^    ^  ^2 

4tt 

Fd    Jl^ 
Yp   =    0    at   X   =    0,    C2    =   — ^ 

Bit 

Fd    £^  ^ 

Yp   =   -4—    (1    -    cos   ^)  (2.9) 


At   X   =    i/2 


Fd   £^ 
■^F  2 
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From  Equation  (2.6) 

K,       d  2  ^  2  ^ 

d    £"*  - 

+   4     (1    -   cos    -J-)  ] 

327T    • 


At   X   =    i/^ 


4  4 

K^  £  d  K^d    I 

v„  =  3141!  tx  ^  -]  -  Tri- 

16tt    EI 


o  F         -^w 

4  4  4 

2  K    £  K    £  K    £    e 

4tt    EI  Idtt    EI 

96F£2         K^^        24K    £4 
384EI   d      =    r    ^        +   — —  +   ^)    d«+   K^    e£ 

TT  2  IT 

4 

.      _   ^if^^ (2.11) 

d      =   J 

°  4       1  24  96F£ 

384EI    -    K      £^     (i.  +   ±|.)     -    ^^Y 

TT  IT 

When  the  effect  of  weight  is  considered: 

4      4 
K;^e£  ±   K^£ (2.12) 

°    384EI  -  0.74638K,£'^  -  96  (-)  F 

1         tt' 

where  K^  =  pAg,  numerically  equal  to  weight  per  unit  length, 


81 


LIST  OF  REFERENCES 


1.  Vanderplaats,  G.  N. ,  "COPES  -  A  FORTRAN  Control  Program 
for  Engineering  Synthesis." 

2.  Vanderplaats,  G.  N. ,  "CONMIN  -  A  FORTRAN  Program  for 
Constrained  Function  Minimization,  User's  Manual," 
NASA  TMX-62,  282,  August  1973. 

3.  Timoshenko,  S.  P.,  and  Gore,  J.  M. ,  "Theory  of  Elastic 
Stability,"  2nd  Ed.,  McGraw-Hill,  1961. 

4.  Timoshenko,  S.  P.,  "Strength  of  Materials,  Part  II," 
3rd  Ed.,  Van  Nostrand  &  Co.,  1955. 

5.  Shigley,  J.  E.,  "Mechanical  Engineering  Design,"  3rd  Ed., 
McGraw-Hill,  1977. 

6.  Tsai,  S.  W. ,  "Composite  Materials  Workbook,"  Technical 
Report  AFML-TR-78-33 ,  Air  Force  Materials  Laboratory, 
March,  1978. 

7.  Rosen,  B.  W. ,  "Stiffness  of  Fibre  Composite  Materials," 
Composites,  January,  197  3. 

8.  Harris,  Bryan,  "The  Strength  of  Fibre  Composites," 
Composites,  July,  1972. 

9.  Barker,  Richard  M. ,  Lin,  Fu  Tien,  and  Dana,  Jon  R. ,  "Three 
Dimensional  Finite  Element  Analysis  of  Laminated  Composites," 
Computers  and  Structures,  Vol.  2,  Pergamon  Press,  1972. 

10.  Puppo,  A.  H. ,  and  Evans,  H.  A.,  "Strength  of  Anisotropic 
Materials  Under  Combined  Stresses,"  Paper  71-368,  12th 
Structures,  Structural  Dynamics  and  Materials  Conference, 
Anaheim,  CA,  April  19-21,  1971. 

11.  Wells,  A.  A.,  "Composite  Materials  and  the  Designer," 
Composites,  May  1972. 

12.  Forest,  J.  D. ,  Hertz,  J.,  and  Christian,  J.,  "Advanced 
Composite  Application  to  Spacecraft  and  Missile  Components," 
SAMPE,  Vol.  3,  No.  2,  January,  1972. 

13.  Vanderplaats,  G.  N. ,  "COMAND  -  FORTRAN  Program  for  Simplified 
Composite  Analysis  and  Design,  User's  Manual,"  NASA  TM-X-73, 
104,  February,  1976. 

82 


BIBLIOGRAPHY 

1.  Loewy,  Robert  G. ,  and  Piarulli,  Vincent  J.,  "Dynamics  of 
Rotating  Shafts,"  SVM-4 ,  The  Shock  and  Vibrations  Monograph 
Series,  Naval  Research  Laboratory,  1969. 

2.  Myklestad,  N.  0.,  "Vibration  Analysis,  McGraw-Hill,  1956. 

3.  Seely,  F.  B. ,  and  Smith,  J.  0.,  "Advanced  Mechanics  of 
Materials,"  2nd  Ed.,  Wiley,  1959. 

4.  Thompson,  W.  T. ,  "Theory  of  Vibration  with  Applications," 
Prentice  Hall,  1972. 

5.  Timoshenko,  S.  P.,  "Strength  of  Materials,  Part  I,"  3rd 
Ed.,  Van  Nostrand  Se  Co.,  1955. 

6.  Vanderplaats,  G.  N. ,  "Numerical  Optimization  -  A  Practical 
Design  Tool,"  Proceedings  ASCE,  7th  Conference  of  Electronic 
Computation,"  St.  Louis,  MO,  August  6-8,  1979. 


83 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Department  Chairman,  Code  69  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

4.  Professor  Garret  N.  Vanderplaats ,  Code  69Vd      10 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Professor  David  Salinas,  Code  69Zc  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

6.  Commandant  (G-PTE-1/72)  2 
United  States  Coast  Guard 

Washington,  DC   20590 

7.  LT  Virgilio  S.  Merced,  USCG  1 
Commander   (mmt) 

Eighth  Coast  Guard  District 
New  Orleans,  Louisiana   70130 


84 


The! 
M49^ 
c.l 


Thesis      ^^^.     189483 

M498   Merced 

c.l        Application  of  numer- 
ical optimization  to 
drive  shaft  design. 


2G  JU^  41 


o  1 
?  1 


I 


Thesi  s 

MW   Merced 

c.l      ^  Application  of  numer- 
ical optimization  to 
drive  shaft  design. 


189483 


thesM498 

Application  of  numerical  optimization  to 


3  2768  001  88605  4 

DUDLEY  KNOX  LIBRARY 


